Quantum Shuttle Phenomena 
in a Nanoelectromechanical Single-Electron Transistor 
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An analytical analysis of quantum shuttle phenomena in a nanoelectromechanical single-electron 
transistor has been performed in the realistic case, when the electron tunnelling length is much 
greater than the amplitude of the zero point oscillations of the central island. It is shown that when 
the dissipation is below a certain threshold value, the vibrational ground state of the central island 
is unstable. The steady-state into which this instability develops is studied. It is found that if the 
electric field £ between the leads is much greater than a characteristic value £ q , the quasiclassical 
£\j ' shuttle picture is recovered, while if £ -C £ q a new quantum regime of shuttle vibrations occurs. 

We show that in the latter regime small quantum fluctuations result in large (i.e. finite in the limit 
h — > 0) shuttle vibrations. 
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i-C , The field of nanoelectromechanics has grown rapidly during the last few years Q, 0, • In particular ,a nanoelec- 
tromechanical sing le-electron transistor (NEM-SET) has been attracting a lot of theoretical |j, la la Lil la 111 llCl 1 1 ll . 
HHQHIIallllllHIiil and experimental [U H3 attention. A NEM-SET is a single electron transistor 
(SET) where the position of a central island (a small metal particle or a single molecule) is not rigidly fixed but can 
oscillate under the influence of an elastic potential. In ]AJ it was shown that in the regime where the island motion can 
be treated classically and the electron tunneling can be described by the Pauli master equation, a new phenomenon 
- a so-called shuttle instability occurs. When a large enough bias voltage is applied between the leads, the island 
oscillates with an increasing amplitude until it reaches a stable limit cycle where it oscillates with some constant 
amplitude. The key issue in |j] was that as the island moves along the classical trajectory its charge q(t) correlates 
with its velocity x(t) in such a way that the time average q(t)x(t) is positive. This results in accumulation of energy 
in the vibrational degree of freedom and in the development of the shuttle instability (see the review [2^j |) 

Further miniaturization of the NEM-SET device brings up quantum effects. In a nanometer-size metal particle, 
fSJ , the electron energy level spacing is about 10 K and the discreteness of the electron energy spectrum can no longer be 
neglected even for temperatures of a few kelvin. In this case the characteristic de Broglie wave length associated with 
the island can still be much shorter than the length scale of the spatial variations of the " mechanical" potential. If so, 
the motion of the island can be treated classically. Shuttle phenomena in this regime have been studied theoretically 
in 0- However, the classical analysis of the shuttle instability (performed in 0) is limited to displacements that 
exceed the amplitude x$ = \Jhj (Mu>) of the zero point oscillations of the island (M is the mass of the island and ui 
its vibration frequency). This quantum limitation raises the question whether or not a threshold value exists for the 
displacement in order for a shuttle instability to develop. To answer this question a quantum theory of the shuttle 
instability must be developed. Moreover, the quantization of the island motion might also effect the steady-state 
regime that develops. 

Different aspects of the NEM-SET in the regime of quantized mechanical motion of the island have already been 
studied H, El El El El El Hj E| • However, no shuttle instability was found because either the coupling between 
electron tunnelling and mechanical vibrations was ignored [a] or strong dephasing in the mechanical dynamics was 
expected El El El El El El- In this Letter we will study the quantum dynamics of a NEM-SET for arbitrary 
dissipation rates. A study complementary to ours has recently been carried out by Novotny et al. (compare |3| and 
|l5l|). However, the numerical analysis reported in Ref. Ilfil was done only for the case when a relatively small number 
of excited vibrational states are involved. This is the case only if the amplitude of zero point oscillations xo is of the 
order of the electronic tunnelling length A* and if the dissipation is large enough. Here we present a complementary 
analytical study valid under the more realistic condition, xq/\+ = A -1 <C 1. 

We will formulate the problem at hand in terms of the dimcnsionlcss displacement x = X/xq and momentum 
p = xoP/h of the island. If we measure all lengths in units of xo and all energies in units of Tito, the Hamiltonian of 
the system reads 

H = ^eakCL^ciak + [e - dx] c ] c + H osc + ^T Q (x) a f ak c + h.c. + H B + H B -osc, (1) 

where a^ fe creates an electron with momentum k in the corresponding lead, a — L, R is the lead index, c 1 creates 
an electron on the single energy level in the island, d = e£ / (Muj 2 xo) is the shift in the equilibrium position of the 
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oscillator due to the electric field £ between the leads, H osc = [p 2 + x 2 ] /2 is the free oscillator Hamiltonian and 
Tl,r(x) — Tl,r{0) exp[=fx/A]. We assume that the electrons in each electrode are non-interacting with a constant 
density of states T> a and that all relevant energies are small compared to the level spacing in the central island which 
for typical systems under consideration exceeds 100 meV. In this case only one single level in the island is relevant 
to the problem. The term Hb describes a heat bath and the last term Hb-osc relates to the coupling between the 
oscillator and the bath 24]. We assume that this coupling is linear in x and treat it in weak-coupling limit. For 
simplicity we will consider only zero temperature case. 

The time-evolution of the system is governed by the Liouville-von Neumann equation for the total density operator. 
After projecting out the leads and the thermal bath we obtain an equation of motion (EOM) for the reduced density 
operator p of the vibrational degree of freedom and the electronic state in the island. Under conditions of large bias 
(eV > hu,e ), the EOM for p becomes Markovian (for details see 0, EE EE IH ) : 



d t p 



[H v + (e - dx)Jc, p] + ttV l (tT l J pdT L - {fgee*, p}) + nV R (2T * R cpJf R - {f^e, p}) + £ 



(2) 



where £ 7 p = — ^ [x, {p, p}] — \ [x, [x, p]], {•, •} denotes the anticommutator, 7 <C 1 is a dissipation rate and time is 
measured in units of uo^ 1 . It follows from Eq. © that the time-evolution of the electronic off-diagonal elements of 
the reduced density operator is decoupled from the evolution of the diagonal elements. After shifting the origin of the 
cc-axis to the point x = d/2 and introducing r a (a;) = 2ttVqT \(x + d/2) we get the system of EOMs for the diagonal 
elements poo =< 1 p > and pn =< l|p|l >, where |1 >= c)\Q >: 



dtpoo ■- 
dtpn 



Hose + Poo 



- {r L (x),p 00 } + \/ t r( x ) PiiV t r( x ) + A Poo , 



H, 



{r_R(x),pn} + y/T L (x) pooV^Ux) + £ 7 pn 



(3) 
(4) 



In what follows we will deal with the evolution of the density operator p + = poo+Pii: which determines the expectation 
values of the observables in vibrational space, and p_ = poo — P11, which describes the shuttling of electrons. 

The problem under consideration can be solved analytically in the limit of a weak electromechanical coupling, 
d/X = e£/(Mu; 2 X i ,) <C 1. This limit has already been used in the classical description of shuttle phenomena |j,|7j], 
where it assured that an electric force acting on the charged island is much weaker than a typical elastic force. For 
simplicity, we will also assume that the tunnelling coupling is symmetric, T^O) = I\r(0) = T/2 and T 1. 

To study the vibrational dynamics near the ground state we use the small parameter A -1 1 to linearize the 
problem with respect to the displacement x. The linearized system of equations, which describes the time-evolution 
of the expectation value ofthe displacement x (t) = (x) and the momentum p(t) = (p), of the island ((•} = Tr{p+(t)»}) 
has the following form 



P, P 



-TP- 



x n_ 

2 



2r_ 

— Tn_ -I — —x . 
A 



(5) 



where n_ = 1 — 2Trpn. An analysis of Eq. (jSJ) shows that an initial deviation from the equilibrium point grows 
exponentially in time with rate constant a — (jthr — 7)/2 if 7 < Tthr — Td/A. Therefore, when the dissipation is 
below the threshold value jthr, the vibrational ground state becomes unstable. 

The exponential increase of the displacement drives the system into the nonlinear regime of the vibrating dynamics, 
where the system may reach a stable stationary state. In order to study this regime we will use Wigner function 
analysis suggested in ^3- The Wigner distribution function (WDF) corresponding to the density operator p± is 
defined by 



W ± (x,p) = ^J dee"** (W||p±k-~> • (6) 
After rescaling the displacement, X = x/X, and momentum, P = p/X, we obtain the following EOMs for the WDFs 



d t W+= Xd P -Pd x +L 1 W+ + L 2 W-, 
d t W- = Xd P - Pd x +t t - r+j W- + \l 2 + r J W+ , 



(7) 
(8) 
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where T± = T R (X) ± T L (X) and 
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(9) 
(10) 



It is convenient to study the steady-state solution in polar coordinates, X = A sin <p, P — A cos p. In these coordinates, 
the steady-state solution is determined by the system of equations 



r+-Li 



'<p-In 

W- = 



W+ = L 2 W- 



l 2 + r_ 



(ii) 

(12) 



with the periodic boundary conditions W±(A, p + 2ir) = W±(A, ip). After eliminating W- from the system of Eqs. 111(1 
and lfH?|) . we get a closed equation for W+ 



0_ : 



L W+ = . 



(13) 



where L = L± + L-z\A — GqL±] 1 Gq\T- + L^] and Go = [dtp + T+] 1 is defined on the space of functions which are 
2-7r-periodic in the variable ip. It is convenient to define a projector V which maps a 2-7r-periodic function f(p) to its 



mean: Vf(ip) = ^ f (ip) dtp / {2tz) and a projector Q = 1 — V . We use these projectors to decompose W + into two 
parts: Vy+(A, <p) = W+(A) + W+(A, ip), where W+ = VW+ and W+ = QW+. By inserting this decomposition into 
Eq. (|13[) and acting on this equation from the left with V and Q, respectively, we obtain two coupled equations for 
W+ and W+: 



PL W+ + W+ =0, 
d v - Ql] W+ = QLW+ , 



(14) 
(15) 



Formally solving Eq. (|15fl for W+ and substituting the result into Eq. I|14(l gives a closed equation for W+(A), 

VL[l-g QL}- l W+{A)=Q, (16) 

where go = c?" 1 acts in the space of 27r-periodic functions with zero mean. One can see from Eq. (|15fl that W+ is of 
lower order in the small parameters d/\, A -2 and 7 than W+. Therefore, in the leading order approximation we can 
write W+(A, ip) w W+(A). If we write the LHS of Eq. (|16|) in terms of A and ip and expand it to second order in the 
parameters d/X, A -2 and 7, we get 



A- 1 d A A[f(A) + D(A)d A }W+(A) = 0. 



where 



f(A) 



D ^4^ + 4A^ 



A 



a (A) 



2A 4 



ai(A) 



d_ 

2A 



1 2 



/3 2 +7^/?3 >0, 



(17) 

(18) 
(19) 



a = -A^VcosipG, ai ee A-^cos^T. <9 p G, G ee G r_ and & = /3 k (A,T). 

We will see later that the functions cto(A) and ceo(A) determine the behavior of W+. The positive function ao(A) 
(see Fig. ^| behaves as T(l + A 2 /2) in the vicinity of A = and for large A it decreases as ln[2A/r]/(7rA 2 ). The 
function a\{A) (see Fig. |2J) is positive and grows as 4A 2 r 3 /9 for small A. For large A it is negative and goes to zero 
as -2/(ttA 2 ). 

Solving Eq. JT7J| gives 



W+(A) = Z^exp. 



/(^) 
I>(A) 



(20) 
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FIG. 2: Function ai(A) for r/(fiw) = KT 1 , 1(T 2 , 10" 3 



where is a normalization constant determined by the condition 2ir J °° dA A W+ (A) = 1 . 

The steady-state solution W + is localized in the phase space around points of maxima of W + . From Eq. f^Hjl one 
can see that W+ has maxima at points Am, where /(Am) = and /'(Am) > 0. In the vicinity of these points, W+ 
is bell-shaped and can be approximated by a Gaussian distribution function with variance a 2 = D(Am) / f (Am)- 
Expanding Eq. (|18(l around A = 0, we find that / ~ A (7 — Jthr)/2 as A — > 0. Thus, W+ always has an extremum at 
A = 0: a maximum if 7 > 7^ and a minimum if 7 < 7t/ lr . This reflects the fact that the vibrational ground state is 
unstable when the dissipation is below the threshold value (the shuttle instability discussed above). 

The global behavior of W+ depends on the electric field £. We have found two different regimes: an electric field 
driven "classical" regime, where £ £ q = C(T)h /(eMA*) and a "quantum" regime, when the electric field is weak, 
£ <C £ q . The dimensionless C(T) ps max[ai(A)]/max[ao(A)] depends only weakly on T (for T — 1CP 3 -7- 10 _1 , it is of 
the order of 10" 2 ). 

In the classical regime, W+ has a maximum at finite A — A c i, if the dissipation is sufficiently weak, 7 < 70 = 
max[ao(A)] d/X > 'jthr- The width of the WDF around A c i is of the order of max{d/A, \~ 3 d~ 1 } 1, which allows 
for a classical interpretation of that regime. The value of A c i corresponds to the stable limit cycle amplitude of the 
classical shuttle oscillations obtained in U This amplitude increases as the value of the dissipation decreases and 
since ag > 0, no stable state with finite A c \ is possible without external dissipation. 

In the quantum regime, the structure of W+ is determined by the quantum fluctuations of the island energy driven 
by inelastic tunneling processes. In this case, the maximum of W+ appears at finite A = A q , when the dissipation 
is below the critical value 71 = max[ai(A)] Xq/(2X+ ). In contrast to the classical regime, A q at low dissipation is 
determined by the zero of the function a.\ (A) (see Fig. |5J| and is still of the order of one even when the dissipation 
is zero. Despite the fact that the amplitude of the shuttle oscillations corresponding to the maximum of WDF is 
much greater then the amplitude of the zero point oscillations, the underlying steady state can not be interpreted as 
classical because the width of VF+(A) around A q is no longer small compared to A q . 

The Wigner function W+ describes the state of the vibrational degree of freedom, while W- relates to the correlations 
between the charge state of the island and the state of the oscillator. It follows from Eq. Ill2[l that W- (A, if) k, 
G(A,(p)W+(A). The WDF for the charged island is given by W\\ = [W + -W-}/2 ss N(A, ip)W + (A), where N = GqT l 
is the occupation of the classical shuttle with an amplitude A and frequ ency ui Q ■ The WDF for the charged island 
exhibits qualitatively the same behavior as was observed numerically in . 

The steady-state current through the system is given by I = eTr[Fi(a;)poo] = f / / dXdPTL(X)[W + + W_] (see 
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[3). Using that W+ w VF+ and W_ « GW+, we get 

/•OO 

7 « (7 C/ (A)} , (.) = 2tt / d J 4 J 4W + ( J 4)[«] , (21) 

Jo 

where /^(A) = ePr^A^GgFi^Y) is the time-averaged current through the classical shuttle 0. Therefore, the 
current / is given by the classical expression I c i{A) averaged over the probability distribution (•) of the amplitude 
A. The function I c i{A) grows monotonically from the tunneling current eT/4 at A = to the shuttle current 
I shuttle = ew/(27r), which is reached already at the amplitudes A s } lut tie — 1. In the classical regime the distribution 
is narrow and I « I c i{A c i), where A c i is the stable limit cycle amplitude of the classical shuttle oscillations. Thus, in 
the classical regime the steady state current is the same as in the classical shuttle case [2j. In the quantum regime, 
A q grows (as the dissipation decreases) only as far as the zero Aq of the function ai(A), but since I c i(Aq) ks I shuttle, 
one can see that the small quantum fluctuations result in a large shuttle current even at E = (as was observed 
numerically in H3). This amplification is another manifestation of the mechanical instability in a non-equilibrium 
NEM-SET. 

In conclusion, we have studied quantum shuttle phenomena in the NEM-SET in the realistic limit, when the electron 
tunneling length is much greater then the amplitude of the zero point oscillations of the island. It is shown that when 
the dissipation is sufficiently low, the vibrational ground state of the central island is unstable. This shuttle instability 
develops into the steady state corresponding to pronounced shuttle vibrations. For large electrical fields between the 
leads this steady-state regime can be interpreted in classical terms. At low field a new quantum regime has been 
found. 
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